Mg NMR study of the MgB2 superconductor 



M. Mali0 and J. Roos, A. Shengelaya and H. Keller 
Physik-Institut, Universitdt Zurich, CH-8057 Zurich, Switzerland 

K. Conder 

Laboratory for Neutron Scattering, ETH Ziirich and PSI Villigen, CH-5232 Villigen PSI, Switzerland 

^^Mg NMR spectra and nuclear spin-lattice relaxation time, Ti, have been measured in polycrys- 
talline MgB2 with a superconducting transition temperature Tc = 39 K in zero magnetic field. From 
the first order and second order quadrupole perturbed NMR spectrum a quadrupole coupling fre- 
quency UQ = 222.0(1.5) kHz is obtained. TiT = 1090(50) sK and Knight shift = 242(4) ppm are 
temperature independent in the normal conducting phase. The '^^Mg Korringa ratio equals to 0.95 
which is very close to the ideal value of unity for s-electrons. The comparison of the experimental 
lyg, TiT, and Kc with the corresponding values obtained by LDA calculations shows an excellent 
agreement for all three quantities. 

PACS numbers: 74.70.Ad, 74.25.Jb, 74.25. Nf, 76.60.Cq, 76.60.-k 



The recent discovery of superconductivity in MgB2 
with remarkably high Tc = 39 K has attracted much 
attention. In particular the observation of a sizeable 
boron Q isotope effect strongly suggests that this simple 
layered intermetallic compound belongs to the conven- 
tional family of phonon mediated BCS superconductors. 
The relevant electron-phonon coupling constant is pro- 
portional to the density of states (DOS) at the Fermi 
level. So it is important to have experimental data of 
this quantity in MgB2. 

Nuclear magnetic resonance (NMR) in metals probes 
the DOS at the Fermi level. The measured quantities, 
the spin lattice relaxation rate, l/Ti, and the Knight 
shift, K, are related to the electron spin susceptibility 
(x(q, uj)) of electrons close to the Fermi level, specifically 
l/(riT) cx lim^^oEqImx(q,t^)/'^, and K a Rex(0,0). 
In case the atomic site symmetry in the crystal structure 
is less than cubic and the atom has a nuclear quadrupole 
moment as it applies for Mg and B in MgB2, quadrupolar 
disturbed NMR delivers in addition valuable information 
on the electric field gradient (EFG) at the specific nuclear 
site. Unlike 1/Ti and K the EFG is determined by the 
distribution of all charges. Up to date only measurements 
of l/Ti, K, and of the quadrupole coupling to the EFG, 
expressed by quadrupole coupling frequency vq, of ^^B 
were performed [|[ 0, |[ While there is reasonable 
agreement among i^g's and relaxation rates, there 
is a considerable controversy concerning the "B Knight 
shifts. While some authors |^ report in the normal con- 
ducting state a small temperature independent isotropic 
shift of 175 ppm, others |^] report a smaller and even 
temperature dependent shift of approximately 60 ppm, 
which they attribute to the Fermi-contact interaction. 
Finally there is also a report |^ of a negative ^^B shift of 
mere 5 ppm attributed to the core polarization. These 
large discrepancies most likely come from the choice of 
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different materials as references for the Knight shift, as 
well as from difficulties to measure the small ^^B Knight 
shift of a broad and in addition quadrupolar ly shifted ^^B 
central line powder spectrum. 

As far as we know, there are no experimental data 
concerning the NMR quantities at the Mg site in MgB2 , 
although theoretical predictions based on ab initio lo- 
cal density approximation (LDA) calculations exist for 
Mg Knight shift, 1/Ti and EFG [@, ||. To some degree 
the lack of experimental Mg NMR data is understand- 
able considering the fact that the only NMR active ^'^Mg 
isotope has a small magnetic moment and low natural 
abundance. Consequently the ^^Mg NMR signals are 
weak even in high magnetic fields. Their detection often 
demands prohibitively long signal accumulation times. 
Despite this handicap we were able to measure NMR of 
the naturally abundant ^^Mg in MgB2. Here we report 
the temperature dependence of Knight shift, spin-lattice 
relaxation time, and quadrupole coupling frequency of 
^^Mg in MgB2. We compare the experimental results 
with the theoretical predictions and we find for all three 
Mg NMR quantities an excellent agreement. Since all cal- 
culations were done by state of the art LDA methods the 
good agreement between calculated and measured NMR 
quantities confirms the LDA as a good approximation for 
MgBa. 

The MgB2 sample was prepared using stoichiometric 
amounts of magnesium and boron (99% and 99.99%, 
respectively) in a form of powder. Both components 
were thoroughly mixed and pressed into pellets. These 
were placed in a tantalum crucible equipped with a non- 
vacuum tight cover. The crucible was then sealed under 
vacuum in a quartz ampoule. The sample was synthe- 
sized during heating at 600, 800 and 900 °C for one hour 
at each temperature. X-ray diffraction (XRD) has re- 
vealed only small amount of MgO as an impurity phase. 
DC magnetization measurements (Fig. |^) in a magnetic 
field of 1 mT (ZFC) yield a transition temperature of 39 
K. 

Before we proceed let us first present a few NMR rel- 
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FIG. 1: Temperature dependence of DC magnetization of 
MgB2 powder in a magnetic field of 1 mT (ZFC). 



evant characteristics of the ^^Mg isotope. The isotope's 
natural abundance is 10.0% and its nucleus has spin 5/2, 
a gyromagnetic ratio 7 = 1.637 x lO^rad s~^T~^, and a 
quadrupole moment Q = 0.1994 x 10"28j^2 g 2by^^ 
NMR measurements were performed on a polycrystalline 
MgB2 powder sample with a pulsed Fourier-transform 
spectrometer at an external magnetic field Bq = 9.047 T. 
Fig. H presents the ^^Mg central line powder spectrum 
gained by Fourier transform of the free induction decay 
signal induced by a 7/Lts 7r/2-pulse. The spectrum shows 
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FIG. 2: ^^Mg central line of MgB2 powder spectrum mea- 
sured at T = 294 K in a magnetic field Bo = 9.047 T. Insert; 
^^Mg first satellites' singularities of MgB2 powder spectrum 
measured at T = 294 K in a magnetic field Bo = 9.047 T. 

the typical powder pattern of the central +1/2 —1/2 



transition line of a spin 5/2 system whose ±1/2 levels 
are disturbed in second order by quadrupole coupling to 
an axially symmetric EFG. There are two peaks in the 
spectrum with an indication of a washed out step left 
and close to the minimum between the peaks. The three 
features from left to right (see Fig. H) are the extrema sin- 
gularities that occur at three powder grain orientations 
having either 9 = 41.8°, 0°, or 90°. Here 6 represents the 
angle between the largest principal component, Vzz, of 
the EFG tensor and the magnetic field direction. Due 
to the site symmetry the EFG tensor at Mg and B sites 
is axially symmetric along the c-axis. Therefore its Vzz 
points into the c-axis direction. From the distance be- 
tween the two peaks in the spectrum and by help of the 
second order quadrupole effect expression |10(] 



Ai/ = 25[/(/-l- 1) - 3/4]- 



we can calculate i/q, defined as 

SeQVzz ^eQVz: 
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Having Aiy ~ 2.85(15) kHz and Larmor frequency 
i^L = 23.570870 MHz, as determined by ^^Mg NMR from 
MgCl2 in water solution, one gets vq — 220(6) kHz. In 
addition we were able to see the singularities of the first 
satellites in the satellite powder spectrum presented in 
the insert of Fig. ||. From the separation of the two first 
satellite singularities, that are positioned symmetrically 
with respect to the central line , one gets immediately 
i^Q. For an axially symmetric EFG vq is simply equal to 
this separation. The gained i/q = 222.0(1.5) kHz agrees 
very well with the above result. The second result, how- 
ever, is always more reliable especially in case when the 
additional Knight shift of the spectrum is anisotropic. 
Since for 9 = 0° there is no quadrupole shift of the cen- 
tral line the 6* = 0° step in the spectrum can be used to 
determine the Knight shift, Kc, of the crystallites whose 
c-axis is parallel to the magnetic field. From the shift of 
the step with respect to i^l we get then Kc = 242(4) ppm. 
Further we noticed that the ratio of the distances of the 
9 = 41.8° and 90° singularities with respect to the step 
is somewhat less than the value 16/9 theoretically ex- 
pected for an isotropic magnetic shift of the central line 
powder spectrum. This allows the conclusion that the 
Knight shift has to be slightly anisotropic. The Knight 
shift, Kab, for crystallites whose c-axis is perpendicular 
to Bq must be a bit larger than Kc- We estimate the 
difference between Kab and Kc to be about -1-10 ppm. 
In the normal conducting phase from 294 K down to 19 
K we observe that the Knight shift (see Fig. |^) and i/q 
= 222.0(1.5) kHz remain in error bar limits constant. 
Such a behavior was also observed for ^^B Knight shift 
ill and J/Q ||, 1). At 19 K the shape of the Mg central 
line spectrum changes drastically . The peak of the 90° 
singularity starts to diminish and disappears quickly by 
cooling below 19 K. In Fig. ^ we exhibit two central 
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FIG. 3: Temperature dependence of the ^^Mg Knight shift, 
Kc, in normal conducting MgB2 powder grains with c || Bq. 
The shift is with respect to the '^^Mg resonance in a water 
solution of MgCl2 at room temperature. The dashed line 
represents the weighted average of Kc- The arrow in the figure 
marks the transition temperature of MgB2 powder grains with 
c ±Bo = 9.047 T . 



line spectra one measured at 6 K and the other one at 
294 K that have their intensities scaled to equal height 
for easier comparison. The peak of the 90° singularity 
seen in the high temperature spectrum obviously is miss- 
ing in the 6 K spectrum. We explain this by the ap- 
pearance of superconductivity in the grains with their 
c-axis close to the direction perpendicular to the exter- 
nal magnetic field. Nevertheless, a substantial part of 
the spectrum coming from crystallites with a smaller 9 
remains unchanged down to 6 K. This indicates that in 
a magnetic field of 9 T at 6 K a substantial part of the 
powder sample still remains normal conducting. In view 
of the large Hc2 anisotropy quoted in the literature (2 - 
9) [0 this comes not as a surprise. However, this has 
consequences for the determination of NMR quantities 
and their temperature dependence in the superconduct- 
ing phase. Necessarily one has to take into account the 
Hc2 anisotropy and avoid to average over the NMR sig- 
nals coming from crystallites with different orientations 
in the magnetic field. For reliable NMR measurements 
in the superconducting state large enough single crystals 
or aligned crystallites of MgB2 are most desirable. 

The measurements of the NMR quantity Ti were done 
by the method of selective inversion of the central line 
and a subsequent monitoring of the nuclear magnetiza- 
tion recovery M{t) at variable delay times t. The re- 
laxation rate defined as 1/Ti = 2W, where W represents 
the spectral density of the fluctuating internal magnetic 
fields, was extracted by fitting the data to the recovery 



FIG. 4: ■^^Mg central lines of MgB2 powder spectra measured 
at T = 6 K (empty circles) and at T = 294 K (filled circles) 
in a magnetic field Bq = 9.047 T. The intensities of the two 
spectra are scaled to equal height for easier comparison . 
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obtained from the solution of the master equation for 
a spin 5/2 system in case of magnetic relaxation and 
selective excitation of the central line . The constant 
r < 2, also a fit parameter, depends on the inversion ratio 
of the central line. We measured Ti only in the normal 
conducting phase at three temperatures 30, 80 and 294 
K. The results are presented in Table |. As expected for a 
metal the product TiT remains constant. The weighted 
average of the three TiT values is TiT = 1090(50) sK. 



TABLE I: Temperature dependence of 
normal conducting MgB2 



^Mg Ti and TiT in 



T(K) 


Ti(s) 


TiT(sK) 


30 


35(1.6) 


1050(50) 


80 


13.5(6) 


1080(50) 


294 


4.35(34) 


1280(100) 



What is the dominant mechanism of magnetic relax- 
ation and Knight shift of Mg and B in MgB2? In most 
metals it is the Fermi-contact interaction of the nucleus 
with the s-electrons at the Fermi level. However, in case 
of MgB2 the states near the Fermi level are mainly boron 
p-electron states with a very small contribution of the 
s-electrons. Pavarini et al. m making ah initio LDA cal- 
culations of 1/Ti and K at "B and ^^Mg sites in MgB2 
show that the dominant relaxation mechanism at the ^^B 
nucleus is the interaction with the electronic orbital mo- 
ment. For the ^^Mg nucleus, however, they predict that 
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the dominant relaxation mechanism is the Fermi-contact 
interaction, which also dominates the Mg Knight shift. 
They get for ^SMg TiT = 1000 sK, Kc = 256 ppm, and 
Kab — 271 ppm which is in excellent agreement with our 
experimental values TiT = 1090(50) sK, = 242(4) 
ppm, and Kab ~ 252(6) ppm. Forming the Korringa ra- 
tio we get experimentally {K'^TiT)/S « 0.95, where 
S = {-iehnf{hlSirfkB = 7.0323 x 10-5 sK for ^SMg, 
with 7e and 7„ the gyromagnetic ratios for electron and 
nucleus, respectively. The observed experimental Ko- 
rringa ratio of 0.95 is very close to the ideal value of 
unity for s-electrons. This result is another confirmation 
that the Fermi-contact interaction is indeed the dominant 
mechanism responsible for relaxation and Knight shift at 
the Mg site. 

As next we compare the experimentally determined 
major principal component, Vzz, of the EFG tensor at the 
Mg site with Vzz calculated ah initio within the density 
functional theory by Tsvyashchenko et al. |^ using the 



full-potential linearized augmented plane wave method. 
From the experimentally determined ^^Mg vq and by use 
of the ^^Mg quadrupole moment Q we extract the abso- 
lute value \ Vzz\ = 3.07(5) x lO^" V/m^ . The quadrupole 
coupling to the EFG and therefore vq do not depend on 
the sign of Vzz- The gained experimental \Vzz\ has then 
to be compared to the absolute value of the calculated 
Vzz = -3.2 X lO^o V/m^ g. We note that the agreement 
of both absolute values is excellent. Tsvashchenko et al. 
g also calculated Vzz at the boron site where they get 
Vzz{B) = 18.5 X 10^° V/m^ which again is in good agree- 
ment with the experimental value Vzz (B) — 17.0(1) x 10^" 
V/m^ ||]. The ability to reproduce the experimental 
NMR quantities at both Mg and B sites by LDA calcula- 
tions certainly strengthens the confidence into the LDA 
approach to calculate other MgB2 material parameters. 

The authors acknowledge useful discussions with M. 
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